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ELECTROTHERMAL OSCILLATIONS AND THE QUASI-LINEAR THEORY OF 
ELECTRON ENTHALPY FLUCTUATIONS IN MAGNETOHYDRODYNAMIC 
GENERATORS AND MAGNETOPLASMADYNAMIC ARC THRUSTERS 

by J. Marlin Smith 
Lewis Research Center 

SUMMARY 

Electrothermal oscillations arising due to fluctuations in electron density and tem- 
perature coupled through a generalized Ohm's Law are studied with regard to magneto- 
hydrodynamic (MHD) power generator and magnetoplasmadynamic (MPD) arc thruster 
applications. The dispersion relation is derived. In the limit of infinite ionization rate 
the dispersion relation reduces to that case which has been extensively analyzed for 
MHD generators. In the MPD thruster configuration this limit results in both a radial 
and a rotational oscillation. For the opposite limit of zero ionization rate, or frozen 
flow, it is shown that an unstable mode can exist which, while in general not critical to 
MHD generator operation, can be responsible for the rotational oscillation observed in 
MPD thrusters. The nonlinear effect of this fluctuation upon the electron properties is 
studied in the quasi-linear limit. Equations are derived for the amplitude of the fluc- 
tuation and its effect upon Ohm's Law and the electron temperature equation. Conditions 
under which a steady state is reached in the presence of the fluctuation are determined. 

The consequences of a purely rotational fluctuation were also investigated. Based 
on the assumption of zero gas temperature in the onset region, which was previously 
used to reduce the existing experimental data for comparison with linear theory, the 
quasi-linear theory predicts unstable growth of fluctuations. Therefore, for the zero 
gas temperature case a "steady state" fluctuation of finite amplitude is not predicted by 
quasi-linear theory. 

For nonzero values of the gas temperature the quasi-linear theory leads to four 
regions in the (Tq /(T g )^ -against- /3 q plane. These regions correspond to absolute in- 
stability (both linear and nonlinear terms are growth terms) , absolute damping (both 
linear and nonlinear terms are damping) , stability at finite amplitude (linear term is 
growth term while nonlinear term is damping), and instability at finite amplitude (linear 
term is damping while nonlinear term is a growth term) . 



INTRODUCTION 


Electrothermal oscillations can arise in gases due to fluctuations in the electron 
temperature. These temperature fluctuations can be amplified through increased ohmic 
heating arising as a result of increased ionization and/or decreased collision frequency. 
These changes in frequency occur as a consequence of the initial electron temperature 
perturbation. The dispersion relation depicting this phenomenon is derived in the first 
section of this report by considering the time -dependent electron density and temperature 
equations coupled through a generalized Ohm's Law. The neutral particles and ions are 
assumed to be incapable of responding at the frequencies of oscillation of interest since 
they are much heavier and much less mobile than the electrons. These heavy particles 
only form a time -independent, uniform background gas in which the electrons fluctuate. 

The dispersion relation is then examined in two limits. In the first limit, the ioni- 
zation rate is taken to be infinitely large compared to the frequency of the disturbance. 
This limiting solution has been widely investigated with regard to MHD power gener- 
ators (refs. 1 to 4). In the geometry of the MPD arc thruster this mode produces both a 
radial and a rotational oscillation. Since this radial component has not been observed 
experimentally in high-pressure MPD arc thrusters, and due to the extensive investiga- 
tion of this mode in the MHD generator case, it is not considered further in this report. 

In the second limit, the ionization is assumed to be frozen. In this limit it is shown 
that the system can become unstable for electron collision frequencies which decrease 
sufficiently rapidly with electron temperature. For the MHD generator this mode does 
not appear to be important since operating conditions and working fluids of interest are 
such that the instability criterion is seldom satisfied. In the MPD arc thruster the gen- 
eral operating range is such that the collision frequency is Coulomb dominated. The 
electron temperature dependence of the collision frequency is such that this mode is un- 
stable and can occur with a purely rotational component, as is observed experimentally. 

In a recent paper (ref. 5) this rotational component was studied in order to make a 
comparison between theoretical and experimental determinations of the rotational fre- 
quency as a function of electric current, magnetic field strength, mass flow rate, and 
atomic weight of the working fluid. Reasonable agreement was obtained even through 
many of the values were far from the stability threshold, where linear results would be 
expected to apply . 

Therefore, in the last section of this report the nonlinear effects of the electron 
enthalpy mode, that is, electrothermal oscillation in the limit of frozen ionization, are 
considered on the basis of a quasi -linear theory. In this theory the governing equations 
are expanded as in the linear case; however, nonlinear terms representing first-order 
corrections are retained. Conditions under which a steady state is reached in the pre- 
sence of the fluctuation are determined. 
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ASSUMPTIONS AND GOVERNING EQUATIONS 


Electrothermal instabilities can arise due to fluctuations in the electron tempera- 
ture. These fluctuations increase the ionization and/or decrease the collision fre- 
quency, resulting in increased electrical conductivity and hence increased ohmic heating 
of the plasma. Under circumstances where this increased ohmic heating occurs, the 
temperature fluctuations are amplified and instability results. 


Assumptions 

In those regimes where electrothermal disturbances are predominant, the following 
assumptions, simplifications, and restrictions are applied: 

(1) The analysis is restricted to sufficiently short wavelengths such that the steady- 
state properties do not vary appreciably over this distance. 

(2) Only fluctuations in space and time of the electron density and temperature are 
considered. This is due to the fact that the relatively lighter and more mobile electrons 
respond to disturbances of greatly different frequencies than do the heavier and less mo- 
bile neutral atoms and ions. 

(3) Only the propagation of magnetohydrodynamic disturbances is considered. This 
restriction implies neglecting the displacement current density in Maxwell’s equations, 
ignoring induced magnetic fields relative to the applied constant magnetic field, and 
assuming the plasma to be qua si -neutral. 

(4) Ion and neutral particle flow velocities are equal; that is, ion slip is ignored. 
Furthermore, all heavy particle flow velocities are assumed equal to the gas flow veloc- 
ity v Q . 

(5) Only propagation in planes perpendicular to the applied magnetic field is con- 
sidered. 

(6) Ion and neutral particle temperatures are assumed to be equal and equal to the 
gas temperature Tq. 

For these limitations, the problem is completely determined in terms of the fluc- 
tuations of electron density and temperature and their coupling through Ohm's Law and 
Maxwell's equations. 


Maxwell's Equations 

From Maxwell's equations, under the assumption of quasi -charge neutrality and 
neglecting the induced magnetic field, the following equations result: 
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V X E = 0 

(1) 


V . j = 0 

(2) 


All symbols are defined in the appendix. 


Electron Continuity Equation 


In the region of interest, the ionization is assumed to be dominated by electron- 
neutral atom ionizing collisions and by three-body recombination. The electron density 
equation then takes the form (see eq. 5. 188 and discussion in chapter 6 of ref. 6.) 




+ X 0 


Vn e = Wi " n e*r 


(3) 


This equation depends upon the steady -state gas flow velocity Vq, rather than upon 
the electron flow velocity v . This dependence arises as a result of assumption 4 and 
the assumption of quasi-charge neutrality. The current density is then given by 
j_= n g e(vQ - v ) and V • j_= 0. Therefore, 


V • 


n,v = V 
e— e 




n e— 0 = lo 


Vn 


(4) 


since the gas flow velocity is constant to within the limits of assumption 1 . 

The electron number density equation is further simplified by assuming that the 
static solution is that of Saha equilibrium so that we are only concerned with small fluc- 
tuations from equilibrium. Under this condition the principle of detail balance (ref. 7) 

can be invoked to relate the ionization coefficient v. to the recombination coefficient v 

i r 

through the equilibrium constant K. Then 

= K> r (5) 


where 


K = 



->3/2 


-T./T 

. i e 


equilibrium 


( 6 ) 
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and T. is the temperature equivalent of the ionization potential (see eq. 6. 83 ref. 6). 
The electron number density equation then has the form 


_9_ 

at 


n 


+ ^0 



( 7 ) 


Generalized Ohm's Law 

The generalized Ohm’s Law is (see ref. 6, eq. 5.168 neglecting ion slip) 



m v 
e 


1 x b 0 


+ -5— VP 

mu e 
e 


( 8 ) 


where E* = E + v q x B q is the electric field in the frame of reference moving with gas 
flow velocity Vq in the applied magnetic field B q. The magnetic field induced by the 
current j has been assumed negligible. The time-dependent and electron inertia terms 
are neglected since the current density relaxes in a mean collision time which is several 
orders of magnitude smaller than the period of the fluctuations of interest. 

The electrical conductivity a is given by 


or 


n e e 


mv 

e 


(9) 


where the total electron momentum collision frequency is 

v=v eK + v ei 

The species collision frequencies are allowed the following general temperature de- 
pendence: 


v ea ^ea n of^e 


ear 


( 11 ) 


where a = i or N, and A and a are arbitrary constants chosen to fit average 

Of W u 

collision frequency data over the temperature range of interest. 
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Electron Energy Equation 


The electron energy equation is 


-k(n - 
2 V e 9t 


T e + n e— 0 


VT 


\ k kT 

7 2 e 


III O ± 

= _2 _® jm — k(T - T) + j • E* 
m a 6 2 


kT. + - kT 
{ 1 2 e 



+ lo 



( 12 ) 


where the VT contribution to the heat conduction term and radiation losses are ne- 
glected. These terms contribute wavelength-dependent damping terms to the dispersion 
relation relative to the wavelength -independent elastic collision damping. The analysis 
is therefore restricted to wavelengths of sufficient magnitude so as to make these terms 
ignorable. The question of the relative importance of these terms has been extensively 
investigated for MUD -generator -type plasmas in references 8 and 9. Equation (12) is 
equation 5.196 of reference 6. , However, q' e in equation 5. 196, as defined by equation 
5.198, is in error. Equation 5.198 is the expression for q g , not q’ e , as can be seen 
from equation 5. 151 and the derivation thereof. Then q’ is related to q by equation 
4. 175 so that the correct value is, neglecting terms the order of the electron-to-heavy- 
particle mass ratio, 


q' = - -P V = -K* 

3-e — e g e — e =e 


VT. 


This is just the portion of the heat conduction vector which by prior argument we ig- 
nored. Therefore, in equation 5.196 of reference 6, V • q* = 0; and our equation (12) 
follows directly. 


UNEAR STEADY-STATE EQUATION 

Under the previously discussed assumptions, the zeroth order (or steady state) so- 
lution is depicted by a region of constant fluid dynamic and electromagnetic properties. 
Then equations (7), (8), and (12) reduce, respectively, to 
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(13) 


2 

- 2 ®- = K(T 
n N0 


eO 


) =K 


0 


lo “ ^0-0 


E* - 


"Vo 


Io x £ 0 


(14) 


i 0 - r o = 2-^o n eof« T eO- T 0> < 15 > 

m a 2 

where the subscript 0 refers to the zeroth order (or steady state) solution. 

The steady state is therefore described by an electron density in Saha equilibrium at 
an electron temperature given by equating the Joule heating to the elastic collisional 
energy losses. 


LINEAR STABILITY ANALYSIS 

The essence of this analysis is to determine whether an infinitesimal perturbation 
to the steady-state solution will initially grow or decay. To this end we therefore con- 
sider solutions of the governing equations of the form 


f = f Q+ f’(r,t) 


(16) 


where fg is the steady -state solution which is constant in space and time and 

f'=real[?e 1( i-£. +a,t> ] (17) 

is an infinitesimal perturbation to this solution such that the products of the f’’s can be 
neglected. The f’s are independent of space and time. We furthermore consider only 
propagation in the plane perpendicular to Bq, which is taken to be in the y -direction, 
so that 


«x- l r l z ) = <J x. < 18 > 

By substituting functions of the preceding form into the governing equations, ne- 
glecting perturbations in the gas dynamic quantities n , v, T and utilizing the steady- 

a — 

state solutions, the equations for the perturbations reduce to the following: 
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col <N 


Maxwell’s equations: 


]x|* =0 

(19) 

l 

a; 

• j_=o 

(20) 

Electron continuity equation: 

/V 

n 

L T. \ T n 


n e0 

( 3 + 1 | e - 2 e 
\ 2 T e0/ T e0 ne0 

(21) 


where 10 * = co + l_ • v q is the frequency in a coordinate system moving with the fluid. 
Generalized Ohm’s law: 


T = a 0 E* + *E*o - ^ I X ^0 + 7 lo 0 + ^ p e 0 1 n T 

m e^0 *0 e^O m e^0 \ eO T e0 


( 22 ) 


where 


n 


/v 

V 


*0 n e0 v 0 


(23) 


= a 


^eNO 3 v eiO\ T e . ^eiO n e 


eN 


+ 


*0 2 V 0/ T e 0 y 0 n e0 


(24) 


and the classical Coulomb potential has been assumed. 


Electron energy equation: 


2 e '^0 n e0 T e0 " T ' 


- ‘(lo • S’ + 1 ' E* 0 ) - (|kT e0 + (25) 
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Dispersion Relation 


The dispersion relation can be obtained formally from the condition for the exist- 
ence of a solution to the homogeneous linear equations (19) to (25), namely, that the de- 
terminant of coefficients be zero. However, since this method cannot be applied to the 
quasi -linear analysis, it is instructive to obtain the dispersion relation by direct sub- 
stitution. This approach is taken in the quasi -linear analysis. 

We start by observing that the perturbed ohmic heating term in equation (25) can, 
with the use of equation (8) , be written as 



The variable E is eliminated from j in equation (22) by forming the vector product 

/V — 

l_x (l_x j) which, by use of equations (19) and (20), results in 



Elimination of E q by using equations (14) and (23) yields 



which upon substitution into equation (26) gives 
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Substituting equation (29) into (25) and using (23) and (24) to eliminate a/o-g and v/vq 
in terms of n e /n eQ and T e /T eQ yields, after some rearranging, 



When equations (14) and (15) are used to eliminate (1/CTq) jg and l_ • ^ j Q x Bq^, the 
preceding equation reduces to 



( 31 ) 
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The dispersion relation is then obtained by using equation (21) to eliminate V*eO from 
equation (31) so that 



Infinite Ionization Rate Limit 


We consider two limiting cases. The first is that in which the ionization rate 
n NO^iO * s muc h greater than the frequency of the oscillation. In this case the electron 
number density is exactly in phase with the electron temperature fluctuation. This is 
seen from equation (21) in the limit O^o^iO 7 ^*^ ~" 00 ’ that is, 


/v 



n e0 



(33) 


This limit is equivalent to assuming that the electrons remain in Saha equilibrium. 
Equation (33) can also be obtained directly from Saha's equation by finding its perturbed 
form. This limit is the one considered by Kerrebrock in reference 1 and has been suc- 
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cessfully used in describing the fluctuations occurring in MHD generators (see refs. 2 
and 3). 

With regard to MPD arc thrusters, this mode was shown in reference 5 to result in 
radial as well as rotational oscillations. Since only rotational disturbances have so far 
been detected in these devices, it is not known how significant a role radial oscillations 
play or even if they exist in these devices. Due to the already extensive investigation of 
this mode, both on the basis of linear analysis and quasi -linear analysis (see review of 
present status in ref. 4), it will not be considered further in this report. 


Frozen Flow Limit 

The other limiting case is the opposite of the case just considered, that is, the limit 
in which the ionization rate * s sma U compared to the frequency of oscillation. 

This limit arises either as a result of the ionization coefficient i^q being small or when 
the gas becomes fully ionized, n^Q — 0. In this case, the electron density is frozen as 
can be seen from equation (21) in the limit that is, 

rsJ 

n_ 

— = 0 (34) 

n e0 


This mode will be called the electron enthalpy mode since only the electron temperature 
fluctuates. Taking the limit (n NO i>, n /a>*) — 0 in the general dispersion relation yields 
for a, = u T - 


0 ) 


* 

r 


— ‘ 2o 

n eO e 


(35) 



The condition that the system be unstable is then 


(i.V\ 

[(i-io) 2 ' 

L "eNO 3 *eiO j 

V T J 

* 2 i o 



(36) 


( 37 ) 
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Since in MHD power generators or MPD arc thrusters T g ^ Tq, the instability criterion 
requires that the total collision frequency be a decreasing function of the electron tem- 
perature (see eqs. (10) and (11) with a oi = 3/2). It is further obvious from the in- 
equality (37) that the unstable disturbance propagates at an angle relative to the steady- 
state current density j Q , and that the growth rate o>* is a maximum for disturbances 
propagating at right angles to j q. 

Physical mechanism for instability of electron enthalpy mode . - The physical mech- 
anism for this instability can be understood by considering a local temperature fluctua- 
tion propagating in the direction as depicted in figure 1 where, for simplicity of argu- 
ment, we shall neglect the Hall effect and the electron pressure gradient in equation (8). 
Consider the region of incremental temperature increase 6T . The total collision fre- 

V/ 

quency depends upon the electron temperature v ~ T* so that 6t> ~ a(6T ). The elec- 
trical conductivity is inversely proportional to the total collision frequency a ~ 1/v so 
that 6a ~ -a(6T e ). In the presence of the electric field, this incremental change in 
conductivity will result in an incremental change in the electric current density. Due to 
the condition V » j = 0 (eq. (2)), this incremental current density is forced to flow per- 
pendicular to l rather than in the direction of the unperturbed current density )q. To 
see how this comes about, we consider the charge continuity equation. The condition 
V • j = 0 arises as a result of the fact that the time for significant space charge buildup 
is inversely proportional to the plasma frequency, which for typical MHD generator or 
MPD arc thruster conditions is many orders of magnitude greater than the 0- to 10 - 
hertz range of interest. Hence, the space charge, for all practical purposes, instan- 
taneously readjusts over the time scale of the wave phenomenon of interest. Thus, over 
the time scale of interest the time rate of change of the charge density term Sq/3t in 
the charge continuity equation (3q/3t) + V • j = 0 can be ignored. Therefore, there is 
an instantaneous buildup of an electric field 6E along l (see eq. (19)) cancelling out 
that component of current density which would otherwise flow across the inhomogenity in 
the _Z direction. To obtain the magnitude of this field, consider the component of change 
in current density along l due to (6a), (6a)(l/Z)£ • E q. The incremental electric field 
in the direction of l needed to cancel this change in current density is then 


(6E) = - ±1(1 • Eq) 

~ "0i ! 


so that the current which actually flows, 6j, is given by 



[lx (6j)] = -(6a) 


— I x (l X E n ) 

o — — — U 
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The resultant change in the local ohmic heating is then given by 


ji 0 • (6E) + (6j) • E 0 = -&2l ± (£ ■ 2 0 )(£ * E q ) - (5a) A- [ft • E Q ) 2 - Z 2 E 2 J 

°0 l 2 l 2L J 


= -a 




l2 io 


‘-Jo)— 

‘"O 7 T e0 


( 38 ) 


where we have used the relations (6a) ~ -a(6T g ) and j q = 0‘ 

For equilibrium to be maintained, the change in ohmic heating must be balanced by 
the change in the collisional energy loss which is given by 


rn , 

2 — - vn a — k(T a - T) 
m 


= (a + 1) 



(39) 


We have used the steady -state relation 

;r^ 2 ^Veof k < T eO- T 0> 

CT 0 m a 2 

and for simplicity have set Tq = 0. The time rate of change in electron temperature is 
then proportioned to the difference between equations (38) and (39) 



which obviously, within the argumentative simplifications made, is fully equivalent to the 
result (eq. (37)) obtained in a formal mathematical way. 
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Having developed this physical understanding of the electron enthalpy mode of insta- 
bility we conclude this section by considering how devices of interest, that is, MHD 
generators or MPD arc thrusters, may be affected by this instability. 

MH D generators . - Two types of MHD generators must be considered. In the 
combustion -type generator, inelastic collisions with the numerous molecular species in 
the combustion gas working fluid make it impossible to raise the electron temperature 
much above the gas temperature . In this case then the factor (1 - (T Q /T e0 )) in the in- 
equality (37) is in general too small for the instability criterion to be satisfied. 

In the nonequilibrium -type generator operating on inert gases seeded with alkali 
metals, the goal is to obtain an elevated electron temperature. An electron temperature 
approximately twice the gas temperature is typical. Therefore, these devices could ex- 
hibit the electron enthalpy mode of instability for values of 



3 ^eio\ 

2 v o) 


-1 


In general, over the temperature range and for neutral gases of interest, the tempera- 
ture characteristic a g ^ is positive or at best only slightly negative; see Maxwell- 
averaged cross sections for electron temperatures near 3000 K presented in refer- 
ence 10. Therefore, it appears that, at least with the proper choice of neutral gas, this 
mode can only become important when the operation is in the upper range of electron 
temperature, where the collision frequency is dominated by electron-ion collisions. 
Furthermore, conditions such as optimal operation tend to require concentration ratios 
of seed (ionizable component) to inert carrier gas (nonionized component) of the order of 
0.1 percent or less (refs. 10 and 11). Since at electron temperatures generally obtained 
(^3000 K) the degree of ionization is of the order of 10 percent or less, the free electron 
concentration is such that collisions with the neutral component are significant if indeed 
not dominant. Therefore, while the electron enthalpy mode could occur in nonequilib- 
rium MHD generators, it appears that its range of importance is small and can be 
avoided without great difficulty . 

MPD arc thrusters . - In order to simplify the discussion of the MPD arc thruster, 
the geometric complexity of the problem is reduced by considering the interelectrode 
gap width to be small compared to the inner electrode radius so that radial variations can 
be neglected. For this situation the geometry can be approximated by two-plane elec- 
trodes parallel to the x-y plane separated by the gap width. The fluid flow velocity and 
applied magnetic field are taken to be in the y -directi on. The applied electric field is in 
the z-direction; and the solution is periodic in the x-direction with period 2nr, where r 
is the mean radius of the fluid flow gap (see fig. 2). For this geometry the steady -state 
Ohm's Law (eq. (14)) has the components 
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^Ox ^O^Oz 


(40) 


where 


•*0y 


= 0 


*0z ~ 


1 + 


% 


?0 = 


eB, 


*Vo 


(41) 

(42) 


(43) 


is the Hall parameter. Substituting equations (40) to (42) into the inequality (37), the in- 
stability criterion becomes 


1 + 2 



3 v _eio\ 
2 *0 / 


< 0 


(44) 


Furthermore, the operating regime of the MPD arc thruster is generally such that the 
electron temperature is much greater than the gas temperature. Also, the degree of 
ionization of the working fluid is generally in the range of 10 percent or better, so that 
the collision processes are dominated by electron-ion collision (refs. 12 and 13), that 
is, (v ei0 /v 0 ) K 1. Therefore, the instability condition for the thruster can be written 
approximately as 


1 



< 0 


(45) 


Obviously, if the disturbance can propagate in any direction relative to E q, the mode is 
always unstable. The maximum growth rate (magnitude of destabilizing or negative term 
is maximum) occurs at the angle given by 


tan 



^1 

^0 


(46) 


where 8 is the acute angle between the direction of propagation ^ and the applied 
electric field IS q. 
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If the mode is restricted by some means to a purely rotational characteristic, that 
is, Z z = 0, then the instability condition reduces to 

|3 q — V 2 "" (47) 

This was the case which was investigated in reference 5, where the theory was applied 
to the onset of the electric -current -formation region near the throat of the device, 
where the pressure is sufficient for condition (47) to be satisfied. 


QUASI-LINEAR THEORY 


We now consider the effect of the nonlinear terms upon the growth of the electron 
enthalpy mode. In doing so, we consider perturbed quantities of the form 

f = (f>+T(r,t) (48) 


The quantity f is the oscillator function for which the time- and/or space -averaged 
value defined by the brackets (f > is zero. In the zeroth order, that is, f = 0, the 
quantity (f) satisfies the steady -state equations (13) to (15). The value of f is taken to 
be sufficiently small so that only first-order corrections to the linear equations need be 
considered. Furthermore, since we are considering small departures, only the funda- 
mental mode is considered; that is, mode coupling between harmonics or other funda- 
mental modes is ignored. Then 

7=f’(t)e i -T- (49) 

Before proceeding to the development of the quasi -linear equations it is appropriate to 
consider some of the average quantities which will appear in the theory. First of all, 
since quantities of interest are all real, definition (49) can be written as 

f = f (t)sini_ • x + f„(t)cosi • x - f 0 (t)sin cp + f„(t)cos <p 
s c sc 

The following pertinent averages over all cp are then obvious: 


17 



The quasi -linear equations are now obtained by utilizing the preceding perturbation 
form and retaining the dominate order in nonlinear terms . 

Maxwell's Equations 

Since equations (1) and (2) are linear in the perturbed quantities, the results are of 
the same form as previously obtained, that is, equations (19) and (20). 

Electron Continuity Equation 

The electron enthalpy mode as previously defined (see eq. (34)) requires that 

n e = 0 (51) 

Generalized Ohm's Law 

Substituting quantities of the form (48) into equation (8) , we obtain 

(|> ♦!-(<*> +”5) (<£*)+£*) -(«>+?)(<!> ♦T)x£+(<<o + (52) 

A 

where b is a unit vector in the direction of B ^ and 

eB 0 

J8 (53) 

“e" 

The average current density is obtained by averaging equation (52) 
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(54) 


<|> =<or>(E*) -<|3><2> Xb-^j)xb + i2 



which we note is just the previous steady-state equation (14) with the addition of the non- 
linear correlation terms x b> and iZ_^<r kT g /e^ • 

We shall now determine these correlation terms as functions of the mean square 
value of the amplitude of the electron temperature fluctuation /t 2 \ . To do this, we 
proceed to determine _j and then (E + i£.kT e /e) as functions of a, / 3 , and the average 
quantities (j) and . From this, (]) is determined as a function of a, /?, and 
^E*^ . In this form the expression is general and applies even in the presence of elec- 
tron density fluctuations. The results are then specialized to the electron enthalpy 
mode by deter mining a and /3 in terms of T . 

The fluctuating current density j is determined by subtracting equation (54) from 
equation (52). 



This equation can be simplified by taking the triple cross product l_x X T) , utilizing 
equations (19) and (20), and using the fact that we are only considering propagation in the 
plane perpendicular to Bq. 


l_x (^x 1 ]) = -1 2 T= vl_x ^l_x ^E*^ - x i x (<j> xb) 


(56) 


The fluctuating electric field E* is determined from equation (55) by taking the dot 
product l_ • J and utilizing equations (19) and (20) to obtain 


<cr>l 2 




~ol_ (l_ • (E*^j (<_J> X b ) + £</3> l • (Jx b) 


+ Higher order terms (57) 


The quantity 



+ i£ kT e /ej will only be used in this report to eliminate the term 
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Y*(p* + *£. ^T e /ejy in equation (54), which is a second-order term. Therefore, the 
terms of second degree and higher in equation (57) are ignored. Substituting equa- 
tion (56) into equation (57) to eliminate j gives 


E* +U— j = -Zl(l • + Ifl • j) xbj 


+ i</3) 


• ((e*) xb)+ pi • o) 


(58) 


The average current density ( j) can now be written as a function of a and p by sub- 
stituting J and E* + i£_kT e /e from equations (56) and (58), respectively, into equa- 
tion (54).. After some amount of algebraic manipulation ( j) can be written in the form 


<|> =<<o(e*) -0X1) xb -.Li 

' -r z 



i- <j> 


(59) 


in which only the dominant order of the nonlinear terms is retained. In this form, the 
equation is quite general and applies to arbitrary electron number density and/or tem- 
per atm e fluctuations. 

These results are now specialized to the case of the electron enthalpy mode by de- 
termining a and 'p in terms of T g when n g = 0. Furthermore, since as previously 
discussed, the instability is in general only important in the Coulomb-collision-dominant 
regime, the subsequent work will be restricted to this limit; that is, — Vq. This 
approximation results in a great simplification of the algebra, without which the subse- 
quent equations become unwieldly. With these simplifications 



m e v 



/v 



(T e > +8 




(60) 


where the subscript 0 in this and subsequent equations signifies evaluation at the aver- 
age electron temperature (T g ) . These quantities will be referred to as zeroth order 
quantities. The average value of the conductivity is then 
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and the fluctuating part 

A# . . 

a = a - (a) 



The quantity of importance throughout the subsequent analysis is the ratio 



(63) 


It is furthermore obvious from definitions (53) and (60) that 



(64) 


The average current density as defined by equation (59) can now be written solely as a 

function of the applied electric field (E* ) , the average electron temperature (T ) , and 

/~ 2 \ 2 ” 

the electron temperature fluctuation amplitude function tT\ /(T e ) through the use of 
equations (61), (62), and (64): 
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Therefore ( j) is completely specified once (T g ) and \TgV/(T ) are established from 
the electron heating equation. However, before proceeding to that task, we place the 
fluctuating current density expression (56) in a simpler and a more convenient form for 
later calculations by noting from equation (59) that 

x<E*)j - <0>i_x[i_x(<2> xb) 

and using equation (64) to yield 


^x(L x <|>) = 


T = — l x(* X ( j>) (66) 

- t 2 (a) ~ V- -> 


Electron Energy Equation 


The electron energy equation (12) can be expanded to yield 

f ta eO^(< T e>^e) = - 2 ^Wo( 1 -f 


O /v» o 

rp mO 

+11 e _ 35 1 e , 


2 <T e> 8 <T e > 2 16 <T e > 3 


^ ~ 2 


x -k((T ) - T 0 M 1 + 5 )+— 1-— + — — 

2 V e on <Te > _ T \ <a) i <a) <ff>2 


a + • • -U(i) 2 + 2(j> • T+ 1 2 

<"> J 


(67) 


Since in the work to follow we shall only be interested in the growth or damping of the 
wave, only terms contributing to this effect have been retained in the preceding equa- 
tion. We also have used equation (8) to replace j • E* by (1/a) j 2 . The advantage of 
writing the fluctuating current density in the form given by equation (66) is now apparent 
since in this form the ohmic heating term in equation (67) reduces to a power series 
in a /{ a ) . 

The average electron temperature equation is obtained by averaging equation (67) 
which results in 
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The equation for the fluctuating part is then obtained by subtracting equation (68) from 
equation (67) 


f kn e0 = - 2 ^ Wo | fc (< Te > - T o) 


3 T e 


/ n ; 

T 


e + 15 I T e 


2 <T P ) <T)-T n 8 




<T P > 




(Tp> 2 / 2 <V- T . v _ e 


e' - 0 \(T„> 2 (T > 2 



e ' ~e' \ / 15 35 

+ 


8<T e >-T 0 16 


X 


< T e>' 


■ + . . 


+ ^_ (j) 2 _J_ 




{<*) 


<J 


<* 2 > N 


<<0 


,<<0 2 <<0 2 


■ + . . . 


<*> 




(69) 


The quantity is determined by multiplying the preceding equation by T e (T e ) , 

utilizing equations (50) and (63) , and averaging to obtain 
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Discussion and Results of Quasi-Linear Theory 

The effect of fluctuations on Ohm's Law. - Of particular interest here is the deter- 
mination of the effective values of the transport properties, that is, the electrical con- 
ductivity and Hall parameter in the presence of the fluctuation as given by equation (59). 
The effective electrical conductivity and the effective Hall parameter are 

defined by the equation 

U > s ff eff<£* > " 0eff<|> x 5 ( 71 > 

Reducing equation (59) to this form yields 



1 + tan^ cp 


where cp is defined as the angle between _Z and (_j_) (0 < cp < 2tt) . By using equations 
(61), (63), and (64), the effective conductivity becomes, in terms of the zeroth order 
conductivity <Tq, 
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(73) 


/p2\ 

1 + 3 il* 


8 <T a > 2 


CT eff “ CT 0 


1 + 


gL 


4 l + tan 2 <p (T > 2 


In figure 3, the ratio of the effective to the zeroth order conductivity is plotted as a func- 
tion of the amplitude function for various angles of propagation. The interesting point to 
note is that the effective conductivity is enhanced by the presence of the fluctuation for 
values of |tan cp | greater than y/E but is reduced for values less than VE. 

The effective Hall parameter is determined by comparing the defining equation (71) 
to equation (59) . The effective Hall parameter is then found to be 


</3> + 


tan cp 


Pefif " 


1 + tan (p 


^+ 




2-) )<£> 

,<*> <p>; 


i + 


1 


1 + tan 2 cp 

<*> 2 


- — ^fy</3> 2 

k <o <js>y / 


(74) 


Using equations (61), (63), and (64), the effective Hall parameter becomes, in terms of 
the zeroth order Hall parameter /3 q, 


Peff ~ 



4 1 + tan 2 cp <T e > 2 


(75) 


Equations (73) and (75) show that in the limiting cases tan cp = 0 and |tan cp | = °° . 


ffeff _ a eff 

% or o 


( 76 ) 
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At intermediate values of |tan <p\, /3 ef f / /3q depends upon the value of /3q. Therefore, 
the dividing point between enhancement and degradation of the Hall parameter is not 
simply a function of the angle of propagation as it was for the conductivity ratio. The 
plot of /3 e g//3Q for the limiting cases is therefore identical to the a e ff/° r o 
ure 3, except that the neutral point /3 e ff / /3 q = 1 occurs for values of tan cp given by 



The values of tan cp are plotted in figure 4. The angles of propagation for which 
j3 e ff = /3 q are restricted to Os tan cp ^ and -°° ^ tan <p =£ - The asymptotic 
values tan cp = ± ylf occur in the limit /3q — 00 . At the limiting point, equality (76) again 
holds, as can be seen by taking the limit /3 q — in equation (75) and comparing the re- 
sults to equation (73) . 

The effect of fluctuations for the case of maximum growth rate . - We are interested 
in determining whether a steady state can exist in the presence of an infinitesimal oscil- 
lation, and if so, under what circumstance the oscillation occurs and what its amplitude 
is as a function of experimentally measured parameters. If the system is steady, then 
3/8t(T e ) = 0 and the balance between Joule heating and electron collision losses as 
given by equation (68) takes the form 


m „ 

2 ^ 1 ”eO>'of k (< T e> - T o) = 


. (i-q>) 2 9 fe 2 ) 

z 2 (j ) 2 ' 4 (V 2 

<V \(*e 2 ) 

8 2 < T e> - T 0/(T ) 2 


<<0 


<i> 


(78) 




This equation reduces to the previously obtained linear steady-state solution (eq. (15)) 
in the limit /(T q ) 2 ~0. 

We first consider the amplitude equation (70) for the direction of propagation which 
results in the maximum growth rate of the disturbance. Since in this direction, the 
fluctuation is first ion stable and grows most rapidly, it is assumed that for small dis- 
placements from the stability threshold this is the observed disturbance. This direc- 
tion, as discussed following the development of the instability criterion (eq. (37)), 
occurs for propagation perpendicular to the average current density; that is, l_ • ( j) =0. 
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For this case, the amplitude function is given by 


1 

1 

m e 

2 _ v 

4-3 T °\ 

+ 3 29, 

( 2 

_ V 

'fit 

I 

$t> 

2 

at <Te> 2 

* "o 

m a 

V < T eV 

2 8 

\29 

T e/ 

'<T e > 2 _ 

<V* 


2 

where we have used equation (78) to eliminate (l/(cr ))(j) and have neglected terms 
greater than second order in ^T^/(T e ) . We note that the linear term is just the sta- 
bility condition obtained from the linear analysis and obviously is damping until the in- 
stability point, Tq/(T ) — 2/3, is reached. Beyond this point the linear term causes a 
growth in the amplitude of the oscillation until the nonlinear term, which is damping for 
Tq/(T > is greater than 2/29, offsets the linear growth term. A steady state is 
reached when 



2 - 



3 29 /_2 
2 8 \29 



(80) 


Obviously, in the case in which Tq/(T ) is less than 2/29 the nonlinear term con- 
tributes to the growth, and the thermal perturbation is unconditionally unstable . How- 
ever, it is obvious from equation (80) that as T o/( T e ) approaches 2/29, ^ T g^/( T e ) 2 
approaches infinity, which violates the quasi -linear condition of small amplitude oscilla- 
tions; that is \Tg)/(T e ) « 1. Therefore, this point is not considered further. 

The amplitude function /(T g ) as given by equation (80) can be calculated by 
determining the temperature ratio Tq/(T ) from equation (68). For l_ • ( j ) =0, the 
steady -state condition 8 /at (T ) = 0, and the geometry shown in figure 2, this equa- 
tion reduces to 



\ < T e>/ 


where 
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2 



( 82 ) 


is the ratio of the ion drift energy to its random kinetic energy and E = (E* ) and Bq 
are the magnitude of the applied electron and magnetic fields, respectively. In deriving 
equation (81) from equation (68), we have also used equations (59), (61), and (64). 

Equations (80) and (81) can be solved simultaneously for given values of the Hall 
parameter (/3) and the parameter R to determine the amplitude function /t^\ /(T ) 
(Note that a physically acceptable solution requires that /Tfi /(T e ) > 0.) These values 

are shown in figure 5, where the amplitude function has been plotted as a function of the 
Hall parameter for various values of the parameter R. The first thing to be noted from 
this figure is that for a given value of R, no physically acceptable solutions exist, that 
is, there is no instability, until a critical value of (/3) is reached corresponding to the 
onset of instability. This critical value of (/3) increases with decreasing values of R, 
as shown in figure 6. Also to be noted from figure 5 is that the amplitude function 
reaches an asymptotic value as (/3) the magnitude of which depends upon the pa- 
rameter R. These asymptotic values are plotted as a function of R in figure 7. 

In figure 8, the ratio of the effective values of the electrical conductivity and Hall 
parameter to their steady-state values are plotted as a function of the zeroth order Hall 
parameter /3 q for various values of the parameter R. These ratios are defined by 
equations (73) and (75), respectively, for <p = tt/2; that is, l_ • ( j) = 0. It is seen that 
the effect of the fluctuation is to enhance the conductivity and Hall parameter, as was 
pointed out previously in figure 3 . 

The effect of fluctuations for the case of a purely rotational disturbance . - In the 
case considered in reference 5 where the propagation was taken to be purely rotational, 
the amplitude equation to second order in ^T^ /( T g ) ^ reduces, after considerable 
algebraic manipulation, to 


1 _a 

2 at 




i 

1 + / 3 ? 



<V 


3 29 /2._Jl0_ + 23 ,32 53 

u/ 2 8 \29 (T ) 29 u 29 (T > 


29 


1+ /3q < T e> 2 


(T p >‘ 


(83) 
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Unfortunately, the dependence of equation (83) upon Tq/(T 6 ) and /3 q is such that it is 
not possible to incorporate the above quasi-linear results in the framework of the linea- 
rized theory of reference 5. The difficulty is that, in order to obtain from the available 
experimental data the values of the local gas dynamic properties required in the theory, 
a number of assumptions had to be made in the analysis of reference 5. One of these 
was to take the local gas temperature Tq to be sufficiently small so that it could be 
ignored in the onset region where the frequency of the oscillation was evaluated. If this 
is done here, the nonlinear term of equation (83) becomes 


iWa + *s.£.j_A j_ 

2 8\29 29 29 % t/! 2 


o 

which obviously is positive and hence destabilizing even for the maximum value of 3 n 

2 ^ 

for which the linear term in equation (83) is destabilizing; that is, (Jq = 2. Therefore, 
a more complete steady-state solution is required and/or local gas properties must be 
measured experimentally if this assumption is the cause of the non -steady -state condi- 
tion. 

However, a second interpretation is possible. Since, in reference 5, the rotational 
frequency was calculated at what was assumed to be the upstream onset point, it could 
well be that the amplitude of the oscillation is growing at this point and reaches its qua si- 
steady -state value downstream of the point in question. In this case, the amplitude may 
indeed be small at the point at which the frequency was calculated, so that nonlinear ef- 
fects can be ignored. The effect of nonlinearity then manifests itself downstream in the 
growth of the oscillation. 

Of further interest is the fact that the linear and nonlinear coefficients in equa- 
tion (83), 


c °" 2 " 3 <v' p “ <84) 

C =329/2.. Jo_ + 23 2 . 53 (85) 

1 2 8 \29 <T e > 29 U 29 U <T e ) 29 / 1 + 

respectively, can be positive or negative, depending upon the magnitude of Tq/(T 6 ) 
and /3 q- In figure 9, the curves of C Q = 0 and Cj = 0 are plotted as functions of 

Tq/ <T e > and /3 q. For both coefficients, points (/3 q, Tq/(T 6 )) above their respective 
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curves result in negative values of the coefficient and hence damp the fluctuation. How- 
o 

ever, for points (/3 q, T^(T e ) ) below their respective curves, the coefficients are posi- 
tive so that the terms contribute to the growth of the fluctuation. As a result, the 

p 

Tq/(T ) -against-^ plane is divided into four regions having the following character- 
istics: 

(1) Region I: Both coefficients are damping, and hence the fluctuation is uncondi- 
tionally damped. 

(2) Region H: The linear term is damping; the nonlinear term is growing. In this 
case the usual linear stability theory would predict a damped disturbance. However, 
the growth contribution of the nonlinear term indicates that, for disturbances with initial 
magnitude of the amplitude function ^T^)/(T e ) 2 greater than |Cq/Cj | , the disturb- 
ances grow unstably . Therefore, in this region, fluctuations are damped for 

/t^) /( T g ) 2 less than |Cq/C^ | and grow unstably for /t 2 \ /(T g ) 2 greater than 

l°o /c il- 

(3) Region IH: Both coefficients are positive, and hence any disturbances grows 
unstably . 

(4) Region IV: The linear term is growing; the nonlinear term is damping. Dis- 
turbance grows or damps until a stable equilibrium point is reached with amplitude 
function ^T^ /(T & > 2 equal to | C Q /C 1 | . 

Whether these regions exist under actual experimental conditions only further ex- 
perimentation can answer. Furthermore, the quasi-linear theory is probably not appli- 
cable to those regions in which it predicts the unstable growth of the disturbance due to 
its requirement of infinitesimal disturbances. In these regions, therefore, a fully non- 
linear theory must be developed. 


CONCLUDING REMARKS 

Electrothermal oscillations due to fluctuations in the electron properties of a par- 
tially ionized gas have been considered. Two limiting cases have been analyzed. In the 
first case, the limit of infinite ionization rate is considered. This limit reduces the dis- 
persion relation to the form which has been extensively analyzed for MHD generators 
(refs. 3 and 4). For the MPD arc thruster configuration this limit results in a rota- 
tional as well as a radially pulsating component. The second limiting case is that of 
zero ionization rate, or frozen flow. For the MHD generator this mode does not appear 
to be important since operating conditions and working fluids of interest are or can be 
chosen such that the instability criterion is not satisfied. However, MPD arc thrusters 
operating at high degrees of ionization are unstable in this mode, and the oscillation can 
occur with a purely rotational component as is observed experimentally. 
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A quasi-linear analysis for the frozen flow limit was also derived. The conse- 
quences of a purely rotational fluctuation for the MPD arc thruster configuration were 
considered. Based on the assumption of zero gas temperature in the onset region, 
which was used in reference 5 to reduce the existing experimental data for comparison 
with linear theory, the quasi-linear theory predicts unstable growth of the fluctuation. 
Therefore, as a consequence of this assumption, a ’’steady state” fluctuation of finite 
amplitude is not predicted by quasi-linear theory. The validity of the zero gas temper- 
ature assumption requires further experimental information or a more complete steady- 
state theory . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 2, 1971, 

112 - 02 . 
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APPENDIX - SYMBOLS 


defined by eq. (11) 

magnetic field vector 

unit vector in magnetic field direction 

defined by eq. (84) 

defined by eq. (85) 

electric field vector 

magnitude of charge of an electron 

electric current density vector 

Saha equilibrium constant 

Boltzmann constant 

wave vector 

particle mass 

particle number density 

pressure 

electric charge density 
defined by eq. (82) 
mean radius of fluid flow gap 
position vector 

temperature of neutrals and ions 

temperature of electrons 

temperature equivalent of ionization potential 

time coordinate 

gas flow velocity vector 

Hall parameter 

electrical conductivity 

total electron momentum collision frequency 
electron-neutral momentum collision frequency 
electron-ion momentum collision frequency 



i>. ionization coefficient 

y recombination coefficient 

r 

(p angle between l_ and j 

6 angle defined by eq. (46) 

oj wave frequency 

-o)j imaginary part of co 

<j) r real part of a> 

Subscripts: 

a heavy particles (neutral atoms and ions) 

e electron 

eff effective value 

N neutral atom 

l_ component parallel to the vector l_ 

x,y,z components of x,y,z coordinate system 
0 zeroth order value 

Superscripts: 

' perturbed quantities 

defined by eq. (17) 

* value in coordinate system moving with gas 

Brackets : 

(...) average value of quantity within brackets 


33 


REFERENCES 


1. Kerrebrock, Jack L. : Nonequilibrium Ionization due to Electron Heating. 

I. Theory. AIAA J., vol. 2, no. 6, June 1964, pp. 1072-1080. 

2. Dethlefsen, Rolf; and Kerrebrock, Jack L. : Experimental Investigation of Fluc- 

tuations in a Nonequilibrium MHD Plasma. Seventh Symposium on Engineering 
Aspects of Magnetohydrodynamics. University of Rochester, N. Y. , 1966, 
pp. 117-125. 

3. Nedospasov, A. V.: Striations. Soviet Physics. -Uspekhi, vol. 11, no. 2, Sept.- 

Oct. 1968, pp. 174-187. 

4. Solbes, Albert: Instabilities in Non-Equilibrium M. H.D. Plasmas - A Review. 

AIAA Paper 70-40, Jan. 1970. 

5. Smith, J. Marlin: Electrothermal Instability - An Explanation of the MPD Arc 

Thruster Rotating Spoke Phenomenon. AIAA Paper 69-231, Mar. 1969. 

6. Sutton, George W.; and Sherman, Arthur : Engineering Magnetohydrodynamics. 

McGraw-Hill Book Co. , Inc. , 1965. 

7. Fowler, Ralph H. : Statistical Mechanics. Seconded., Cambridge Univ. Press, 

1966. 

8. Lutz, Michael A. : Radiation and Its Effect on the Nonequilibrium Properties of a 

Seeded Plasma. AIAAJ., vol. 5, no. 8, Aug. 1967, pp. 1416-1423. 

9. Hiramoto, Tatsumi: Critical Condition for the Electrothermal Instability in a 

Partially Ionized Plasma. Phys. Fluids, vol. 13, no. 6, June 1970, pp. 1492- 
1498. 

10. Heighway, JohnE.; and Nichols, Lester D. : Brayton Cycle Magnetohydrodynamic 

Power Generation with Nonequilibrium Conductivity. NASA TN D-2651, 1965. 

11. Smith, J. Marlin; and Shair,F. H. : Theoretical Optimum Seed Concentrations in 

Slightly Ionized Nonequilibrium Plasmas. Phys. Fluids, vol. 9, no. 8, Aug. 
1966, pp. 1615-1617. 

12. Workman, J. B. : Arc Thruster in a Magnetic Annular Discharge. AIAA J., 

vol. 7, no. 3, Mar. 1969, pp. 512-519. 

13. Rosciszewski, Jan: Acceleration Process in an MPD Arc -Jet. Tenth Symposium 

on Engineering Aspects of Magnetohydrodynamics. W. D. Jackson and J. E. 
Klepeis, eds., Avco Everett Research Lab. , 1969, pp. 112-116. 


34 



Figure 1. - Depiction of electron enthalpy fluctuation. 



(a) Assumed geometry, 



(b) Approximate geometry used for analysis. 


Figure 2. - Thruster geometry. 





igure 3. - Ratio of effective to zeroth order electrical conductivity as 
function of amplitude function for various angles of propagation. 


tan 



0o 


Figure 4. - Value of tan<p at which effective Hall parameter equals 
zeroth order Hall parameter (3g as function of p 0 . 







Figure 9. - Damping and growth regions for rotational dis 
turbances. 
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